Strong effect of weak diffusion on scalar turbulence at large scales 
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Passive scalar turbulence forced steadily is characterized by the velocity correlation scale, L, 
injection scale, I, and diffusive scale, r^. The scales are well separated if the diffusivity is small, 
Vd <C l,L, and one normally says that effects of diffusion are confined to smaller scales, r <^ r^. 
However, if the velocity is single scale one finds that a weak dependence of the scalar correlations 
on the molecular diffusivity persists to even larger scales, e.g. I ^ r ^ ra We consider the 
case of L 2> i and report a counter-intuitive result - the emergence of a new range of large scales, 
L ^ r ^ P /rd, where the diffusivity shows a strong effect on scalar correlations. 



Studies of passive scalar advection in a random smooth 
flow were pioneered by Batchelor [1] and Kraichnan Q 
who considered the opposite extremes of almost frozen 
and short-correlated in time random velocity gradients. 
The two approaches were later extended into a unified 
theory describing the statistics of scalar correlations in a 
general smooth flow These theoretical stud- 

ies were originally motivated by interest in explaining 
the so-called viscous-inertial interval of advection at the 
scales smaller than the viscous, Kolmogorov, scale. How- 
ever, the theory, which is nowadays often called Batchelor 
flow theory, also applies to many cases of non-turbulent 
but chaotic smooth flows, e.g. of the type discovered 
recently in polymer solutions j8| . 

Main theoretical efforts in the field were focused on 
the analysis of scalar correlations within the convective 
range, ^ r <C i.e. at the scales smaller than 
the injection scale / but larger than the diffusive scale 
Td [1, M The range of scales above the pumping 

scale, even though very nontrivial with highly intermit- 
tent correlations [ll|, attracted much less attention. In 
this work, we continue to discuss the domain of large 
scales. Complementary to our general interest in un- 
derstanding multi-point correlations in turbulence, this 
study was additionally motivated by our recent interest 
in the condensate regime of 2d turbulence [l2|, where 
small scale vorticity is advected passively by the large 
scale, coherent part of the flow. 

On a superficial level studying correlations of a fluc- 
tuating quantity upscales from the injection scale may 
seem akin to many problems in equilibrium statistical 
mechanics, e.g. of the type considered in the field of crit- 
ical phenomena where one studies fiuctuations of an order 
parameter driven by thermal noise at small scales. How- 
ever, the essential difference here is that our problem is 
off-equilibrium. A particularly important consequence of 
this fact is an intermittent, strongly non-Gaussian statis- 
tics observed for the problem both at the scales smaller 
d, Q and larger [llj than the pumping length I. 

In this work we extend the analysis of Ref. [llj ac- 
counting for the effects of molecular diffusivity which 
were ignored in (llj . A surprising result of our study is 
that diffusion, even though small, dominates correlation 
functions of the scalar at large scales, r ^ P /vd- This re- 



sult is Batchelor fiow specific and it can be explained in 
dynamical, Lagrangian, terms. The coUinear anornaly, 
established in [9| and later discussed in @, [l^, [H, Ej], 
states that Lagrangian particles released along a line in 
a Batchelor fiow stay aligned, unless weak diffusive ef- 
fects are accounted for. The anomaly reveals itself in an 
angular singularity of the passive scalar multi-point cor- 
relation functions observed near the parallel alignment of 
the points [§, [l3| • Translation of the coUinear anomaly 
from the dynamical to statistical language goes as fol- 
lows. If diffusivity is neglected, a blob of freshly injected 
passive scalar is deformed by a smooth flow into a strip 
of the same density. The strip contributes to correlation 
function of the passive scalar provided it covers all the 
points where the correlations are measured. Thus the 
strip should have the length r of the order of the sepa- 
ration between the points and it should also be oriented 
in a way that all the points are covered. Since the flow 
is chaotic and orientation of the stripe is random, the 
probability to cover the points is determined by the an- 
gular size of the stripe. In incompressible flow, the blob 
volume is conserved. The volume can be estimated as 
Z"*, where d is the space dimensionality, thus the cross- 
section of the stripe can be estimated as l"^ /r and the 
angular size of the stripe is ^ jr'^ . This results in the 
cx scaling for the n-tli order correlation function of 
the scalar, if = {6i ■ ■ ■ On), in the coUinear geometry, 
i.e. when the points v\, • • ■ , Vji lie on a straight line, and 
r is the size of the most separated pair of points [ll|- 
Volume preserving stretching of the scalar blob, injected 
at the pumping scale I, should be modifled when the blob 
size in the contracting direction reaches ^ r^, since the 
diffusion blocks further contraction of the blob beyond 

Td- 

Let us concentrate on the 2d case. For a seriously 
stretched stripe, with the spatial extent r > P/rd, 
the stripe is ~ r^-wide in cross-section and therefore 
(Xd ~ Td/f gives an estimate for the angular size of the 
stripe. The temporal dynamics of the stripe is as follows. 
The stripe grows in size (along the stretching direction) 
while the scalar density inside the stripe, estimated as 
^ Td/r fraction of the initial density, decreases. Besides, 
different stripes stretched simultaneously start to over- 
lap along the contracting direction. Due to the random 
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FIG. 1: Figure Caption. Schematic plot illustrating La- 
grangian (temporal) evolution of two blobs of scalar. (A) 
Initial injection. The blobs are of size ~ I separated by the 
distance ~ I. (B) Result of diffusion-less deformation. The 
blobs grow in size along the expanding direction of the flow 
and decrease in size along the contracting direction. Vol- 
umes and initial concentrations of the scalar inside the blobs 
are preserved. This phase terminates when the width of the 
blobs reaches r^. (C) Further deformation keeps the width 
of the blobs ~ r^, while the lengths of the blobs continue to 
increase. Volumes of the blobs grow and the density of the 
scalar inside the blobs decreases. Blobs will eventually over- 
lap. Dashed lines show projected shape of the blobs realized 
if the diffusivity is ignored (naive extension of the stage B). 



character of the passive scalar injection, the sign of the 
density in the overlapping stripes alternates, thus leading 
to destructive interference. This additional effect leads 
to further suppression of the scalar correlations by the 
factor where, N ^ rvd/P, is an estimate for the 

number of stripes (that were initially of size I and sepa- 
rated by the distance ~ I) which contribute to the over- 
lapped conglomerate. Combining the pieces, one derives 
the following scaling for the 2n-th order correlation func- 
tion of the passive scalar measured at r l'^ /vd within 
the collinear geometry: 



(1) 



Eq. ID) is the main result of the paper, that will be con- 
firmed below with a proper rigor. The result shows a 
strong dependence of the high order correlations of the 
scalar at the scales beyond the injection scales on the 
diffusivity, and it should thus be contrasted with a much 
weaker dependence on diffusion observed in the passive 
scalar correlations at the scales smaller than the injection 
scale [9]. 

Even thought Eq. ([T|) is derived in 2d, the qualitative 
result, stating a strong sensitivity to diffusivity of the 
scalar fluctuations at large scales, also extends to 3d (and 
higher dimensions). In general, the simultaneous correla- 
tions are expressed in terms of the Lagrangian evolution 



of a fluid blob, while diffusivity stops contraction of the 
blobs at the diffusion scale, r^, thus making the blob glob- 
ally sensitive to the small scale, diffusion-related physics. 
Obviously, the Sd picture of the phenomenon is more 
evolved due to existence of an additional, third, dimen- 
sion in the blob dynamics that can be either contracting 
or expanding. As a result, the 3d generalization of Eq. ([T|) 
becomes sensitive to the sign of the second Lyapunov ex- 
ponent of the flow. This sensitivity is similar to effects 
discussed in [l^ and [ij] in the contexts of kinematic 
dynamo, and decaying scalar turbulence, respectively. 

The dynamic equation for a passively advected scalar 
held, 0, is 



(2) 



where u{t,r) is the flow velocity fleld, n is the diffu- 
sion coefficient and r) is the pumping term. The 
velocity u and the forcing 4> are assumed to be indepen- 
dent and random functions in space/time with prescribed 
statistics, spatio-temporally homogeneous and spatially 
isotropic. The forcing is correlated at the scale, /. We 
consider Batchelor, that is spatially smooth, flow where 
the velocity fleld is correlated at the scale L, the largest 
scale in the problem. We also assume that the velocity 
fluctuations are sufflciently intense to guarantee that the 
diffusive range, r <C r^, where the effects of advection are 
strongly suppressed by diffusion, is realized at the scales 
smaller than the pumping length, i.e. < I. In Batche- 
lor flow, the velocity difference between points separated 
by a distance much smaller than L is given by the flrst 
term of the Taylor expansion in the inter-point separa- 
tion, Ua{ri) - Ua{r2) ~ <Jaf}{ri-p ~ 'r2-p)- Therefore, the 
velocity derivatives matrix, a, is the only velocity related 
characteristic entering the problem at the scales smaller 
than L. In an incompressible flow, discussed here, the 
velocity gradient matrix is traceless, tr tr = 0. We also 
assume that a, followed in the reference frame of a fluid 
parcel, is finitely correlated in time. 

Representing solution of Eq. ^ in the Lagrangian 
frame (see [1| for derivation details) one arrives at the 
following formal expression for the scalar field 



6)(t,r)= / dt' exp 




Ht',R)- (3) 



Here R = W{t',t)r and W is the ordered exponential 



W{t',t) ^Texplj dr a{T) 



(4) 



Note that in an incompressible flow, deiW{t) ~ 1. The 
argument R{t') of the function in Eq. ^ traces back 
in time the Lagrangian trajectory arriving at the position 
r at the moment of time t. The K-dependent exponential 
on the right-hand side of Eq. ([3]) represents effects of 
diffusion. Therefore, Eq. ([3]) is nothing but a formal 
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way to express the aforementioned quahtative arguments 
concerning the Lagrangian evolution of a passive scalar 
blob. Since (/) is spatially correlated at the scale I, the 
temporal integral on the right-hand side of Eq. ^ is 
formed at i — i' ~ X~^ln{r/l) where A is the principal 
Lyapunov exponent of the flow, defined as the average 
logarithmic rate of Lagrangian trajectories divergence. 
This stretching time diverges as r ^ oo, which allows 
us to consider the forcing as short-correlated in time. 
Then the forcing field is effectively Gaussian and thus 
fully described by the pair correlation function 



X{r-r')= dt m,r)c^{t\r')), 



(5) 



decaying sufficiently fast with increase in r at r > L We 
assume that x('") is a function of |r| only. 

For the short-correlated forcing any correlation func- 
tion of the scalar field can be calculated in two steps. 
First, one averages over times larger than the pumping 
correlation time but smaller than A^^ ln(r/Z). This is for- 
mally equivalent to averaging over the statistics of forc- 
ing for a given realization of the velocity field. Averaging 
over velocity, corresponding to longer times and larger 
spatial scales, follows. This scheme gives the following 
expression for the simultaneous pair correlation function 
of the scalar, K{r) = {e{t,r)9{t,0)), 



lC{r) 



dt I dk (exp(J)xfe) , 



(6) 



J ^ikW{-t,Q)r~2nkI{t)k, (7) 

i{t)^ I dTW{-t,-T)W^{^t,~T), (8) 

Jo 

where Xk is the Fourier transform of x(r) and T indicates 
matrix transposition. The only averaging left to be done 
in Eq. ([6]) is over statistics of a. 

Consider the d = 2 case and introduce an Iwasawa-like 
decomposition for the ordered exponential 



cos </J sm Lp 
— sin cos ip 











1 C 

1 



(9) 



This representation is useful as the three governing fields, 
(yj, Q, C decouples in the asymptotic limit of large time, 
t ^ A^^. Moreover, at the large times, the orientation 
angle p becomes random uniformly distributed over the 
range (0;27r), the distribution function of C, freezes to a 
non-universal stationary shape and the typical C becomes 
a fluctuating 0(1) value, while the probability distribu- 
tion of the finite time Lyapunov exponent, A = g/t, at- 
tains the following self-similar form [ll| 



V{t, X) (xVt exp [-tS'(A)] 



(10) 



Here <S'(A) is the so-called Cramer function which is con- 
cave and achieves its minimum at A = A, and the ^/t 
factor accounts for accurate normalization of the total 



probability to unity. In the limit Xt ^ 1, the main con- 
tribution to the integral ^ originates from r close to t, 
thus leading to 



lit) 



.2e(t) 



o{t) 
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d-\t), 



(11) 



where c is a fluctuating factor of order unity and O is the 
93-dependent part of the decomposition ([9]). (See Q for 
detailed discussion of the c-field statistics.) Averaging 
over homogeneous random orientations (p (reflecting the 
assumed isotropy of the velocity fluctuations) one derives 
the following expression for J from Eq. ([7]) 



J = ir (fciCe^ + k2e~^) - 2crlkle^e ^ 



(12) 



where ~ k/X and ki,k2 are components of the wave 
vector k in the reference frame fixed by the decomposi- 
tion dH). 

A comparison of the two terms in Eq. (I12p suggests 
that the outer scale interval, r 3> splits in two distinct 
sub-intervals: I <^ r <^ l^/rd and r ^ l^/r^^. To describe 
the first interval of relatively small scales one may ignore 
the last term in Eq. (|12p . Then, direct integration of 
Eq. ([6]) results in the diffusionless scaling 



/C(r) cx— d'^x x{x) 



(13) 



already derived in [ll|. However, evaluating the integrals 
in Eq. ^ in the regime where the second term in Eq. (fT2|) 
dominates the first one does not actually change the final 
answer for the pair correlation function ()13p . Indeed, in 
this limit, integration over fci in Eq. ^ is determined by 
the diffusive exponential, that allows simply to replace fci 
by zero in the integrand of Eq. ^ . Integrating the result 
over t one arrives at a factor S{k2)/r, while subsequent 
integrations over k2 and over the domain of small C ^ 
1, leads to the same expression for JC{r), independent of 
the diffusion coefficient. 

We will see now that the cancelation of the r^j de- 
pendence in the pair correlation function is incidental, 
and it does not actually extend to the general case of 
higher order correlation functions. Consider, for exam- 
ple, the fourth order simultaneous correlation function, 
K(^Hri,r2,rs,ri) = {e{t,ri)e{t,r2)6{t,r:,)e{t,ri)) , 
which is decomposed into the following sum, 
C(ri2,r34) + C(r23,ri4) -I- C(ri3,r24), in the case of a 
Gaussian pumping, where rab = fa — ri,- Being interested 
in establishing scaling law for the special case of coUinear 
geometry, one focuses on analysis of C(r,r). Generaliz- 
ing evaluations resulted in Eqs. (|6ll2p . one arrives at the 
following expression valid at r ^ I, 



C{r, r) oc 



dti 



exp 



(fciCe^ 



^26 



-2crl {kl 
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where Qi = g{ti) and g2 = gih)- If Z <C r <C P/rd then 
the diffusive exponent in can be neglected and one 
arrives at the diffusionless expression 



relation functions cancel. In this case the behavior of 
the correlation functions is determined by non-universal 
features of the flow statistics. 



C{r,r)(x dh dt2 (x {rCe^' , re'^' ) x He^^ , re'^^ ) ) 



leading to the scaling C(r, r) cx derived in [TT|. Note 
that the main contribution to the above time integrals 
comes from the region exp(pi) ^ exp(£i2) ^ r/l ^ 1. In 
the r ^ P /rd limit, the diffusive exponential in Eq. ([M]) 
cannot be replaced by unity. On the contrary, it domi- 
nates integration over ki and qi resulting in emergence 
of the X(o,fc2)X(o,i32) term in the integrand. Then, inte- 
grations over ti and t2 decouple from each other and one 
arrives at 



C(r, r) cx 



1 



(15) 



in accordance with the general formula ([T]) . 

The strong dependence of the correlation function on 
the diffusion, observed for the coUinear geometry, does 
not extend to a general off-collinear case, where thus 
the diffusionless consideration of [Tll| applies. These dis- 
tinct coUinear and off-collinear results are asymptoti- 
cally matched in the r^/r-small angular vicinity of the 
coUinear geometry. 

Note also that if the Corrsin integral / d^r x{f) is equal 
to zero then the leading terms and in the cor- 



Summarizing, we have shown in this work that weak 
molecular diffusion does control the large scale correla- 
tions in scalar turbulence steered by the Batchelor incom- 
pressible flow. The main logical points of this work are: 
(a) Correlation of the passive scalar within the coUinear 
geometry are much stronger than in an off coUinear case. 
The angular extent of the coUinear anomaly domain is 
controlled by the fact that a scalar stripe injected and 
stretched by the flow cannot get thinner than rd- (b) The 
effect of diffusivity leads to a faster decay of scalar corre- 
lations with the scale r at the largest scales than in the 
domain of smaller scales, where diffusion is irrelevant, (c) 
Scaling in the diffusion-controlled regime becomes sen- 
sitive to the order of correlation function and (in 3d or 
higher dimensions) on the character of the flow statistics. 
The situation in 3d, even though qualitatively similar in 
the sense of a general importance of effects of diffusion, is 
actually more involved because of an interplay between 
two distinct Lyapunov exponents characterizing the flow. 
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